Reversal of Klein reflection in bilayer graphene 
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Abstract 

Whereas massless Dirac fermions in monolayer graphene exhibit Klein tunneling when passing 
through a potential barrier upon normal incidence, such a barrier totally reflects massive Dirac 
fermions in bilayer graphene due to difference in chirality. We show that, in the presence of magnetic 
barriers, such massive Dirac fermions can have transmission through even at normal incidence. The 
general consequence of this behaviour for multilayer graphene consisting of massless and massive 
modes are mentioned. We also briefly discuss the effect of a bias voltage on such magnetotransport. 
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One of the most remarkable features of electron transport in graphene is that charge 
carriers are chiral in nature ^ and their degree of chirality changes as the number of layers 
are varied [3- 5] . In monolayer graphene, chirality results in Klein tunneling, namely the 
perfect transmission of a normally incident electron of energy E through a potential barrier 
of height V when E < V. In contrast to this, for bilayer graphene chirality of a different de- 
gree leads to total reflection from the same barrier at normal incidence, a phenomenon often 
ca.le d Klein Section Q Q. Tflns, « Stt on gly influences coherent flaiflsfc transnns- 
sion through graphene based heterostructures, a phenomenon which is now experimentally 



accessible [8l-lll| and potent with possibilities for new device applications. 

For monolayer graphene, it has been shown that the combined effect of an inhomogenous 
magnetic field, dubbed as a magnetic barrier, and an electrostatic potential significantly 
alters ballistic transport of charge carriers. It has already been shown that use of such 
inhomogenous magnetic barriers can lead to confinement of chiral electrons as opposed to 



Klein tunneling 



12| . This has subsequently led to a large body of work [l^ • 
In this Letter, we consider the effect of such barriers in the case of bilayer graphene. 
Fabry- Perot like transmission fringes develop due to scattering of chiral charge carriers from 
inhomogeneous magnetic barriers in the presence of voltage. By analyzing these fringes, we 
find that, in contrast to the above discussion on monolayer graphene, it is quite the opposite 
in bilayer graphene. It is seen that it is possible to actually reverse the Klein reflection 
of normally incident electrons; namely, the barrier can cause transmission and significantly 
modifies the observed conductance in bilayer graphene. As we show, this also paves the way 
of a more generalized understanding of the effect of such magnetic barriers on transport of 
chiral quasiparticles in multilayer graphene which shares many features of monolayer and 
bilayer graphene. 

The generalized chiral operator of degree J that describes chirality of charge carriers 
in multiple layer graphene with ABAB • • • type of stacking is h J v — nj • cr, where rij = 
cos(J(f))x + sin( J(j))y, (f) = tan -1 ^ and a is the two-component Pauli matrix. Monolayer 
and bilayer graphene are just two cases of this description with J = 1 and J = 2 respectively. 
The general Hamiltonian for charged quasiparticles of J-th degree chirality is then 

Hj cc \p\ J h J p (1) 

The above Hamiltonian leads to a pseudo Zeeman effect due to a momentum-dependent 
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FIG. 1: E-k diagram and electron trajectory in bilayer graphene with (a) magnetic barrier, and (b) 
barrier + well 

pseudo magnetic field B p = \p\ J n [jjj]. The corresponding eigenstates with definite chirality 
of degree J are given as 
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(2) 



These eigenvalues are proportional to ±\p\ J for the conduction (+) and the valence (— ) 
bands respectively and give both massless and massive Dirac modes for multilayer graphene 
in general |5|. Eq.([2j) also demonstrates how the pseudospin winding number described in 
Ref . fis| varies with the number of layers. 

The upper and lower components of the pseudospinor give the probability amplitude 
of an electron being in sublattice A and B. A pure scalar potential V(r) cannot flip the 
pseudospin, namely (±|]/(r)|=F) = 0, since it couples in the same way with both sublattices. 
For monolayer graphene with J — 1, there is only one massless Dirac mode. Since a scalar 
potential shifts the Dirac point, to preserve pseudospin, an electron (hole) state outside the 
potential barrier needs to be matched with a hole (electron) state inside the barrier and that 
leads to perfect transmission. 

The present work discusses bilayer graphene in Bernal stacking, which corresponds to 
J = 2 in Eq.(pp) and has one massive Dirac mode. Here also, a scalar potential cannot flip 
the pseudospin and shifts the Dirac point. However, now the spectrum is parabolic. Hence, 
a state with p 2 outside the potential barrier will go to a state with — p 2 inside the barrier to 
preserve the pseudospin. This requires that a propagating wave outside the barrier should 
be matched with an evanescent wave inside the barrier leading to total reflection at normal 
incidence. 



3 



Unlike a scalar potential, a vector potential due to a magnetic field B = B{x ) y)z couples 
with the momentum and can flip the pseudospin. In the Landau gauge the vector potential 
is A = Ay{x)y ) one can show that when such a magnetic field is present the pseudospin 



projection angle changes from Jcj) to J(j) r where (\J = tan _1 (^-), where q y = k y + ^ (cf. 
Fig{T](a)). Here, we use k = {k X) k y } with no magnetic barrier and q = {q X) q y } with a 
magnetic vector potential. The eigenstates are then 



\+)b = 



V2 



e 2 

• j ($) 
e 2 
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(3) 



Scattering between one of the eigenstates in Eq.([2j) with B = to one in Eq.([3]) for B ^ 
is now possible as the pseudospin is not necessarily conserved in the presence of a magnetic 
field. Consequently, perfect Klein reflection upon normal incidence will not take place in 
bilayer graphene in presence of such a magnetic barrier. We show this by explicitly calcu- 
lating transmission by a transfer matrix approach Q ■ Chiral charge carriers in bilayer 
graphene obey 
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(4) 



Here, they are described by a 4-component spinor i$!{x,y) = (^r fl ty b ty c ty d ) T and a 4 x 4 
Hamiltonian in the presence of a magnetic barrier and an electrostatic potential, n = 
Vf[Px + i(p y + e A/c)] with Fermi velocity vp — 10 6 m/s. The magnetic field profile is 
B = BQ(x 2 — d 2 )z, Vi and V2 are the potentials at the two different layers, and t is the 
tunnel coupling between the two layers. This Hamiltonian reduces to a manifestly chiral 
symmetric form given in Eq.(pp) in the limit j < 1, where E is the incident electron energy 
m- Vi = {¥^2 correspond to unbiased (biased) graphene bilayer cases. The stationary 
solutions \t = ^(x,y)e - ^ which obey Hip = Eip can be explicitly obtained in the three 
regions x < — | (region I), — | < x < | (region II) and x > | (region III). In regions I (j = 1) 
and III (j = 2) these solutions are propagating as well as evanescent waves. In region II 
where the magnetic field is finite, stationary solutions are given by parabolic cylindrical 
functions. We define £b = an d £b = jrf as units of length and energy respectively 

and introduce dimensionless variables x —¥ , V\ i2 = -7^ and e — In regions I and III, 
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the dispersion relations now become 

[-q{xf + (e + S) 2 ][-q(x) 2 + (e - S) 2 ] = (e" - 5 2 )f 

Here, q(x) 2 = q 2 x + q 2 y , q y (x) = k y + sgn(x)it§-^ e = e - v with v = (v x + v 2 )/2 = ^ 
and S — (vi — v 2 )/2 = Also, is the total magnetic flux through an area Mb in units 
of the flux quantum $ = and £b = |§ is the magnetic length. For an evanescent wave, 
Q 2 X = ~4 < o. 

We now begin with the case of an unbiased graphene bilayer where 5 = and Vi = v 2 = v . 
In regions I and III we set = Vq = 0. The transmission through such a combination 
of magnetic barrier and electrostatic barrier can be written using transfer matrices A4q(x) 
and Mb{x), where the subscript (B) defines the region with magnetic field (£>). These 
matrices are given below explicitly. 



M B (x) = 



D p+ {z) D p -{z) D p+ {-z) D p -(-z) 

elp+Dp+^z) elp'Dp-^z) D p +_ r (-z) e\p~ 'D p -_ r (-z 

etD p +{z) s 2 D p -(z) £%D p +(-z) £ 2 D p -(-z) 

£ 1 £%D p ++ 1 (z) £ 1 £2D p -+ 1 (z) £{£%Dp++ 1 (-z) 8^ D p - z) _ 



Mv{x) = - 



e f e -i>QxX 



e'e~ KxX 



[q x - iq y (x)]e iq * x -\q x + iq v (x)]e iQxX i[n x - q y (x)]e KxX -i[n x + q y (x)]e KxX 



-e'e iq * x 



-e'e~ iq * x 



e'e KxX 



e'e KxX 



(5) 



_-[q x + iq y (x)]e iqxX [q x -iq y (x)]e iqxX i[k x + q y (x)]e KxX -i[k x - q y (x)]e KxX 
Here, z = V2(x+k y ),p ± = ( 7± -l)/2,^ = ^ and 4 = e ¥ - and 7± = e' 2 ± y/l + e f2 t' 2 . 



The current density expression is obtained as j x = 



ip. The transfer matrix 



through any combination of a scalar and vector potential can now be written in terms of 
transfer matrices given in Eq.([5]) for regions with finite and zero magnetic field. This can 
then be used to find the transmission. The ratio of current density in region III and the 
incident current density in region I gives the transmission probability as a function of the 
angle of incidence cj) = tan -1 

In Figl2l the plots in the left column show transmission through a magnetic barrier with 
increasing strength as a function of the incidence angle and the applied voltage in the barrier 
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FIG. 2: Transmission T as a function of incidence angle for different fields, (left) Barrier, (right) 
Barrier-hwell. 

regime. The uppermost figure shows the transmission for £> = 0, which clearly shows the 
region of perfect reflection at and around the normal incidence symmetrically placed between 
two wings of resonant Fabry-Perot fringes. This is a generic feature of transmission through 



such barriers 



a. 



As the strength of the barrier increases, a transmission region develops 
between these two wings due to the effect of magnetic barrier and the resulting transmission 
also becomes asymmetric. Clearly seen is the disappearance of perfect Klein reflection at 
normal incidence. The total angular range of transmission however shrinks in the presence 
of the magnetic barrier since, beyond a certain incident angle, all electron waves suffer total 
internal reflection. This is an important result of our paper since it shows the exact reason 
why an inhomogeneous barrier causes Klein reflection to be reversed to transmission for 
bilayer graphene. 

In the right column in Figj2] is plotted similar transmission through a setup consisting 
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FIG. 3: G for different magnetic barrier strengths, (a) Barrier, (b) Barrier+well. A — 0, E — 
17meV kept constant. 

of two magnetic barriers, equal in magnitude and opposite in direction such that the total 
flux through the region vanishes, as shown in Fig{TJ Now, the incident and the transmitted 
wave vectors are parallel to each other. As a result, the net rotation of the pseudospinor 
due to the inhomogenous field vanishes and Klein reflection at normal incidence is restored. 
The bending of the Fabry- Perot fringes can be attributed to the asymmetry in the angular 
transmission at a given voltage in presence of magnetic barrier + well combination. 

We now study the effect of the above transmission on conductance at very low temper- 
atures and for energies close to the Fermi energy in the linear trans por t regime. To a good 
approximation, the dimensionless conductance G can be written as |13] 

/ZL 
2 #T(£,0)cos0 (6) 
-I 

A comparison between the conductance with and without a magnetic barrier in Figj3] shows 
that the gaps between the conductance maxima and minima get reduced in the presence of 
a magnetic barrier. The barrier reflects electrons incident upon it beyond a critical angle 
and the angular range of transmission shrinks with increasing barrier strength. As a result, 
the absolute value of conductance maxima comes down. For the barrier+well case, the 
conductance is slightly higher. This detailed characterization of ballistic transport through 
magnetic barriers dictated by pseudospin chirality is one of the main results of this work. 

The preceding discussion shows how the effect of localized magnetic field on transport 
of massive Dirac modes encountered in bilayer graphene is very different as compared to 
that of massless Dirac modes in monolayer graphene. Another important observation can 
je made by looking at magnetotransport through magnetic barriers for monolayer graphene 
12| and the present analysis for bilayer graphene. The same barrier which is reflective for 
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FIG. 4: Biased bilayer transmission for various V+, A and magnetic field strengths (colours and 
lines same as in Fig{^). 



the massless mode is transmissive for the massive mode. The general theory of N-layer 
graphene predicts the existence of 1 massless and N — 1 massive modes for TV odd and only 
massive modes for TV even j|J. The magnetic barrier thus can be used to at least partially 
filter out or to selectively allow a mode. 

Finally, we briefly discuss the influence of magnetic barriers on transport through a biased 
bilayer; i.e., when V\ ^ V2 in Eq.(|4j). The effective 2x2 Hamiltonian is 



4(n) 2 



(7) 



H u \) is the Hamiltonian for the unbiased bilayer given in Eq.(pQ) for J = 2. Thus, bias 
introduces an additional z component of the pseudo magnetic field. This will try to rotate 
the pseudospin out of the x — y plane. If A is large and the incident energy lies in this gap, 
the pseudo Zeeman gap will suppress transmission. This is because pseudospin-flipping is 
energetically costly, a phenomenon that has been experimentally observed [181-42 if] . The third 
term in Eq.([7j) does not flip the pseudospin. Combining the first two terms, the Hamiltonian 
can be rewritten as 



(n • a 



'A 2 



h 2 q 2 Vp X 



Here, a are three-component Pauli matrices and n = (sin 6 cos 20', sin 6 sin 20', cos 6) with 



tan# 



h 2 q 2 v 2 F 
tA ' 



The corresponding pseudospinors |+) = (cos |e l( ^\8m^e l ^) T and |— ) 



(sin |e - ^', — cos |e^') T are again eigenstates with definite chirality. Comparing these states 
with the ones given in Eq.([2j), we find that the magnetic barrier gives an in-plane rotation 
while the bias voltage gives an out-of-plane rotation to the pseudospinor. The two orthogonal 
rotations can also be given in different spatial regions and have been recently discussed using 
a Berry phase argument [22] for a pure bias. In this paper, we have considered the case 
where A and B are non-zero in the same region. If the gap is large and the incident energy 
lies inside the gap (FigJH(a)), the transmission as expected is highly suppressed and there is 
no visible effect of the magnetic barrier. Also, we find that, as long as the incident energy 
does not lie in the gap (FigJH(b), (d)) or the gap is too small (FigJH(c)), the effect of the 
barrier on the biased and the unbiased bilayers is similar. 

This work is supported by grant SR/S2/CMP-0024/2009 of DST, Govt, of India. One 
of the authors (NA) acknowledges support from CSIR, Govt, of India. 



[1] A.H. Castro Neto et al, Rev. Mod. Phys. 81, 109 (2009). 

[2] M.I. Katsnelson and K.S. Novoselov, Sol. St. Comm. 143, 3 (2007). 

[3] K. S. Novoselov et a/., Nat. Phys. 2, 177 (2006); 

[4] E. McCann and V.I. Fal'ko, Phys. Rev. Lett. 96, 086805 2006. 

[5] H. Min and A.H. MacDonald, Phys. Rev. B 77, 155416 (2008); H. Min and A.H. MacDonald, 

Prog. Theor. Phys. Suppl. 176, 227 (2008). 
[6] K.S. Novoselov et al, Nat. Phys. 2, 620 (2006). 
[7] S. Park and H.-S. Sim, Phys. Rev. Lett. 103, 196802 (2009). 

[8] A.V. Shytov, M.S. Rudner, L.S. Levitov, Phys. Rev. Lett. 101, 156804 (2008); M. Ramezani 

Masir, P. Vasilopoulos and F.M. Peeters, Phys. Rev. B 82, 115417 (2010). 
[9] A.F. Young and P. Kim, Nat. Phys. 5, 222, (2009). 
[10] N. Stander, B. Huard, and D. Goldhaber-Gordon, Phys. Rev. Lett. 102, 026807 (2009). 
[11] A.F. Young and P. Kim, Ann. Rev. Cond. Matt. Phys. 2 1-20 (2010). 
[12] A. De Martino, L. Dell' Anna, R. Egger, Phys. Rev. Lett. 98, 066802 (2007). 
[13] M. Ramezani Masir, P. Vasilopoulos, A. Matulis and F.M. Peeters, Phys. Rev. B 77, 235443 
(2008). S. Ghosh and M. Sharma, J. Phys. Cond. Matt. 21, 292204 (2009); M. Sharma and S. 
Ghosh, J. Phys. Cond. Matt. 23, 055501 2011; L. Dell'Anna and A. De Martino, Phys. Rev. 



9 



B 80, 155416 (2009); L Dell' Anna and A. De Martino, Phys. Rev. B 79, 045420 (2009) ; Y.P. 

Bliokh, V. Freilikher and F. Nori, Phys. Rev B 81, 075410 (2010) 
[14] A.K. Geim and A.H. MacDonald, Phys. Today 60, 35 (2007) 
[15] C.-H. Park and N. Marzari, arXiv:cond-mat/1105.1159 (2011). 

[16] M. Ramezani Masir, P. Vasilopoulos and F.M. Peeters, Phys. Rev. B 79, 035409 (2009) ; M. 

Ramezani Masir et al., App. Phys. Lett. 93, 242103 (2008). 
[17] J. Nilsson, A.H. Castro Neto, F. Guinea and N.M.R. Peres, Phys. Rev. B 76, 165416 (2007) ; 

M. Barbier, P. Vasilopoulos, F.M. Peeters and JM. Pereira, Phys. Rev. B 79, 155402 (2009). 
[18] T. Ohta et al, Science 313, 951 (2006) 
[19] E. V. Castro et al, Phys. Rev. Lett. 99, 216802 (2007). 
[20] J. B. Oostinga et al, Nat. Mater. 7, 151 (2008). 
[21] Y. Zhang et al, Nature 459, 820 (2009). 
[22] S. Park and H.-S. Sim, arXivxond-mat/1103.3331 (2011). 



10 



T 



^ o 



-2 



